Spectroscopy for cold atom gases in periodically phase-modulated optical lattices 
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The response of cold atom gases to small periodic phase modulation of an optical lattice is 
discussed. For bosonic gases, the energy absorption rate is given, within linear response theory, by 
the imaginary part of the current autocorrelation function. For fermionic gases in a strong lattice 
potential, the same correlation function can be probed via the production rate of double occupancy. 
The phase modulation gives thus direct access to the conductivity of the system, as a function of 
the modulation frequency. We give an example of application in the case of bosonic systems at zero 
temperature and discuss the link between the phase- and amplitude- modulation. 
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Cold atomic systems have proven to be remarkable lab- 
oratories to study several effects of strongly correlated 
systems. In particular the control of parameters, kinetic 
energy in an optical lattice and interaction using a Fes- 
hbach resonance, allows us to potentially use them as 
quantum simulators, with considerable success both for 
pure and disordered systems [l|, |2j. However, in addi- 
tion to realizing the systems, the ability to probe it is 
important. Because of the electrical charge neutrality of 
cold atoms, unlike electron systems, they are insensitive 
to the usual electromagnetic probes. This makes it po- 
tentially difficult to probe correlations in such systems. 
To overcome this issue, several probes have been pro- 
posed besides the standard time of flight (TOF) exper- 
iment such as Bragg spectroscopy [3rl6| to measure the 
dynamic structure factor, radio frequency spectroscopy 
measurement 0, Q to count the number of molecules 
formed by the Feshbach resonance, shot noise measure- 
ment P-KLljl for the density-density correlation function 
or momentum-resolved Raman spectroscopy EDI for 
the single-body spectrum function. 

Among the various spectroscopic probes a particularly 
simple probe consists in ch anging periodically the ampli- 
tude of the optical lattice [LJ, |l5[ . The energy absorbed 
by such a modulation can be estimated from the TOF im- 
age. The correspondingtheory of the energy absorption 
rate (EAR) spectrum [16], [TtJ , was shown to give access 
both to the Mott-insulating (MI) gap and to the kinetic 
energy correlations in the system. Although measuring 
the EAR by the TOF was possible for bosons, a simi- 
lar measure was highly inconvenient for fermions. It was 
proposed [18| that a measurement of the doublon pro- 
duction rate (DPR) in response to the amplitude modu- 
lation would give access to the same information. Such a 
measure was successfully implemented for fermionic sys- 
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The amplitude modulation of the optical 
lattice coupled either to EAR or to DPR is thus a simple 
but powerful and versatile probe. 

In this Letter we propose an alternative probe, based 
on a phase modulation of an optical lattice potential. 
Such a modulation is known to lead to a current 123142a 



or to band narrowing [271 . |28| . Here we use the phase 
modulation in connection with EAR or DPR techniques, 
to analyze the spectrum of the system. We show that 
such a probe gives access to the current autocorrelation 
function and is thus analogous to optical conductivity 
measurements in condensed matter systems, allowing a 
very close comparison at the experimental level between 
the two domains. We illustrate the use of such a probe by 
some examples for bosonic gases and compare with the 
spectrum obtained by the amplitude modulation spec- 
troscopy. 

Let us first describe our proposed probe: The optical 
lattice potential is created by shining laser against a mir- 
ror. If the mirror is stationary, the created D-dimensional 
optical lattice is given as V op (r) = Vq Xm=i c °s 2 {Q p. r 
where Q = (Qi, • • • , Qd) is a wave vector of the optical 
lattice. One can modulate the phase of an optical lattice 
potential by oscillating the mirror as shown in Fig. [T] 
The lattice potential in the laboratory frame is modified 
as V op {r,t) = E^i^ocos 2 ^^ - F M (t))] where F(t) 
represents the oscillation of the phase. It is convenient 
to switch to the comoving frame by the gauge transform 
U(t) = exp(iMF(t) • J/h) where M is a mass of the 
atoms and J the current operator [29| . In the comov- 
ing frame, the lattice becomes a stationary one, V op (r), 
and an additional term, which reflects the inertial force, 
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V op (x,t) = V cos 2 [Q(x - F(t))] 



F(t) 



FIG. 1. (color online). A schematic showing the setup of the 
periodic phase modulation of an optical lattice. The incident 
laser and the reflected one forms the standing wave which 
corresponds to an optical lattice. The lattice potential follows 
mirror oscillation, and consequently the phase is modulated. 
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emerges in the Hamiltonian: 



H(t) =H - MF(t) ■ J, 



(1) 



where H is the Hamiltonian of the interacting system in 
the optical lattice V op (r). Carrying out the further gauge 
transform U'(t) = exp[i/VM fdt' F 2 (t')/2h] where N is 
the total atom number, we can find an expression iden- 
tical to that of a charged particle in an electromagnetic 
field, the atomic mass corresponding to the charge. The 
vector potential is given as A ext = F(t) with V-A oxt = 
and the scalar potential is zero cxt = [53] . The system 
thus behaves as charged particles under an electric field 
E e xt(t) — —F(t). Hereafter we set F(t) — /cos(wi). 

Let us first consider bosonic atom cases. One can 
then measure the EAR by similar techniques than for 
the amplitude modulation 
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The EAR is given 
by the time average of the absorbed energy: R(fj) = 
^ f^ +2n ^dt{H(t)) where (• • • ) denotes the statistical 
average with the Hamiltonian ([T]). The EAR within the 
linear response theory is given by 



^pm(w) — — 
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where ffi^, (w) is the Fourier transform of the retarded 
current correlation function for the Hamiltonian Hq. 
Note that for the EAR, due to H(t) = Muj 2 cos(wt)/ • J, 
one can automatically derive the second-order response 
of R(lu) in terms of within the first order perturbation 
theory. Note, as shown in Eq. (JTJ, that in order to stay 
within linear response a small modulation is necessary. 
In particular one needs <C Qn 1 , thus a modula- 

tion amplitude smaller than a lattice constant. This is 
something difficult but achievable with the current ex- 
perimental technique |30j . For the higher frequency the 
smaller amplitude needed to stay within the linear re- 
sponse and away from the dynamically induced phase 



transition [3l|, |32j. We will confine our analysis in the 
following to such a regime for which the EAR gives di- 
rect access to II^(cj). 

Since the EAR directly gives I$ v (ui), it is immediately 
related to the "optical conductivity" [33| . While the zero 
frequency part of the Drude peak in the conductivity will 
be suppressed in the EAR due to the factor u 3 in Eq. , 
all the other features, at finite frequency, are perfectly 
reproduced. It is thus a particularly useful quantity to 
make comparison with similar phenomena in condensed 
matter systems, or to probe the physics of disordered 
systems, for which transport is the prime probe. 

It is interesting to compare this result to another one 
obtained for the amplitude modulation |16(, in the same 
linear response regime. In the later case the EAR is either 
given by the density correlation function for weak opti- 
cal lattices or by the kinetic energy one for strong optical 



lattices 16|, |22( . The different representation of pertur- 
bation operator comes from the fact that the energy scale 
of the amplitude modulation also goes beyond the chem- 
ical potential as the energy scale of the lattice potential 
V op (r) increases. On the contrary, in the case of the 
phase modulation, the perturbation in Eq. (JT|) does not 
follow the energy scale of the lattice potential. Therefore, 
Eq. ((2|) is independent of the strength of the lattice poten- 
tial. This is a definite advantage of the phase modulation, 
which is always related to the same physical quantity irre- 
spectively of the strength of the optical lattice. Another 
important difference between the phase and amplitude 
modulation comes from their symmetries. Indeed, for 
example, in the case of a ID strong lattice the amplitude 
modulation would correspond to the kinetic energy cor- 
relation function which is given by T cx cos(k)b\.bk 
while the current is given by J cx ^ fc sin(fe)&tfrfc' For 
the parity inversion, J > — J while T is invariant. This 
affects the selection rules. The phase modulation pertur- 
bation probes the transitions from a state to an oppo- 
site parity state. In contrast, the parity is preserved in 
the transition due to the amplitude modulation. Thus, 
both modulations complement each other in the way they 
probe excited states, and lead in general to different re- 
sults. 

For fermions, as for the amplitude modulation, the 
EAR is not a convenient way to probe the conse- 
quences of the modulation. We follow here the same 
approach as in [18] and show that for the phase mod- 
ulation the measurement of the DPR gives essentially 
the same information as the EAR. We assume that Hq 
in Eq. ([I} is described by the Hubbard model, Hq = 
~*H E CT ,(i,j) C L C ^ + UJ2j n 3t n ji- The number of dou- 
bly occupied sites is defined as Njj(t) — (V)/U, and can 
be rewritten as N B (t) = ^ [(#(*)) - (T) + MF ■ (J)], 
where T and V are the kinetic energy and interaction 
terms in the above Hamiltonian. The production rate 
is defined as the time average of Nn(t) for a period: 
P(oj) = £ p+^^dt'Nnit'). A second-order perturba- 
tion expansion in F(t) gives for the productions rate 



-PpmM = R P m(u)/U, 



(3) 



which shows the direct relation between the EAR and 
the DPR for the phase modulation. This shows that DPR 
gives also access to the optical conductivity for these sys- 
tem. 

The results @ and the equivalence of the DPR to the 
EAR ([3]) are thus our main proposal to use the phase 
modulation of the optical lattice to measure the optical 
conductivity of interacting systems in a cold atom con- 
text. 

Let us now examine an example of the phase modu- 
lation technique. For the sake of simplicity we take a 
repulsively interacting ID bosonic atom system at zero 
temperature. The unperturbed Hamiltonian in Eq. ([1} is 
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written as 



Hq = dx 



• 2M x 



» + V op (x) V + ^P 



(4) 



where V op (x) = Vocos 2 Qx is a ID lattice potential, and 
g an interaction parameter. The field ip( x ) an d p( x ) are ) 
respectively, the annihilation and density operators. 

For a shallow lattice potential, Vq <C /z, V op (x) can be 
regarded as a perturbation, and then the Hamiltonian (|4]) 
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can be rewritten via the bosonization 
yfpe~ lB ( xS> andp(x) ~ p—d x (p(x)/n+2pcos{27rpx—2ip(x)], 
where p is the mean density of the system. Retaining 
only the most relevant term generated by the presence of 
V p(x), one can obtain 

H cS = H TL + A dx cos [2(Q - irp)x + 2<p(x)] , (5) 



where Htl = ^ Jdx [K(d x 9(x)) 2 + R-^dMx)) 2 ] is 
the Tomonaga-Luttinger (TL) Hamiltonian. The phase 
fields 9{x) and ip(x) represent, respectively, the phase and 
density fluctuations of bosons. For an arbitrary repulsion 
K runs from oo to unity as the interaction increases, and 
K = 1 and oo correspond to the Tonks gas and the nonin- 
ter acting bosons, respectively. Thus in the boson systems 
the low-energy physics is mainly governed by the cosine 
term in Eq. Furthermore, the current is written as 
J = vK Jdx U(x) where H(x) is the canonically conju- 
gate momentum to <p{x). In the incommensurate case, 
i.e., irp ^ Q, the cosine term in Eq. ([5J vanishes, and the 
effective theory is the TL liquid. For Tip = Q, the model 
([5]) becomes a sine-Gordon model for which the cosine 
term is relevant for K < 2, leading to a MI gap in the 
excitation spectrum 35[. We will consider this last case 
in what follows. 

The conductivity can be calculated using the methods 
in Refs. [H, |36j|. To determine the large frequency be- 
havior of the ID current correlation function ITj(u;), we 
use the memory function method, which gives correctly 
IIj(cj) for large frequency compared to the MI gap. The 
memory function M(uj) = wn^(w)/[IT}(0) - n^(cj)] can 
be approximated as M(u) « [n£(0) - Uf(w)]/wU^(Q) 
where 11^ (uj) is the retarded correlation function of 
F(t) = [H ,J(t)}. From the Hamiltonian ©, F is 
given by F — i2vKX Jdx sin[2(Q — up) + 2ip(x)], and 
M(u) ~ lu 2K ~ 3 is immediately computed. In the gap- 
less (K > 2) case, negligible M(w)/w for small w leads 
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As a result, it is found that the 



EAR spectrum for small ui and K > 2 behaves as 
Rpm(w) oc uj 2K ~ 2 . This is to be compared with the am- 
plitude modulation result [16| Ram{w) oc uj 2K ~ 1 for weak 



lattices. A similar result is obtained in the large-w limit 
for the massive case K < 2. In the gapful case the cosine 
is relevant, and the conductivity, i.e., the phase mod- 
ulation response will be zero below the gap 35]. This 
example thus shows differences between the phase and 



amplitude modulation. This difference has two origins; 
one is the trivial different prefactors of the correlation 
functions (w 3 for the phase modulation and w for the am- 
plitude one). More importantly and as discussed above, 
the main difference comes from the perturbation cou- 
pling to two different operators: namely the current for 
the phase modulation and the density for the amplitude 
modulation in the shallow lattice limit. 

We now consider a strong lattice potential. Then, the 
system is well described by a lattice model: Hq is given by 
the Bose-Hubbard Hamiltonian H c s = T + V where T = 



h.c] and V = % J2 3 n j ( n j ~ !)■ an d 
are, respectively, the creation and number operators 
for a bosonic atom at the jth site, and in and U are the 
hopping parameter and on-site interaction, respectively. 
For an incommensurate filling, the ground state is in the 
gapless superfluid (SF) phase. For filling of one particle 
per site, the SF-MI transition occurs at U/ta = 1-92 [37| . 
In the SF phase, the low-energy physics is governed by 
the TL liquid [35j. In the MI phase, an energy gap opens, 
and the low-energy physics is no longer described by the 
TL liquid. The lowest energy excitation above the gap is 
formed by a pair of atoms at the same site (doublon) and 
an empty site (holon). In the limit tn/U — > 0, the pair 
excitations are Af(ftf— I)-fold degenerate where is the 
number of lattice sites. For finite but small tn/U, the 
degenerate energy levels split, and an energy band whose 
width is about in is formed. This band leads in the phase 
and amplitude modulation spectrums to a peak around 
w w (7/fi as shown in Fig. [2j 

Using degenerate perturbation theory [16(, the EAR 
spectrum can be calculated. Let \d^h r ) be a pair 
state of the doublon and holon at i?th and rth 
site, respectively, which is an exact eigenstate of 
V. We represent T in the Hilbert space spanned 
by \dnh r ). The eigenstate of T is \l,l') = 
W Ea=i EI'i 1 e i(p,R+ar6K ' lf) sm( qi ,r)\d R h R+r ) where 
w Pl = 1 + 2e tp ', vi = 2-KljM and q v = irl' /Af (I = 
1, • • • ,Af and /' = 1, • • ■ ,Af — 1). The corresponding 
eigenenergy is Ei t \i = U — tn\w Pl | cos (ft/. The retarded 
correlation function of O for w > is expressed as 
3ng(w) = -TrhJ2 n \(n\O\0)\ 2 6(huj - E n ) where \n) and 
|0) = nj^l uac )' respectively, denotes an intermediate 
state and the MI ground state. Restricting the inter- 
mediate states onto \l,l'), both the amplitude modula- 
tion and the phase modulation can be computed. In the 
Af — > oo limit they turn out to be identical and can 
be written as 3TIj (aj) 
-Afht H Tr(u}) 



-Af^-n(uj) and 3Ll£(w) 

4 



[16( where tt(cj) = fyl - [(fiw - U)/3t H } 2 , 
as shown in Fig. [2j Therefore, the appropriately scaled 
EAR, i.e., w _2 i?pM(w), is identical to Aam(w). In this 
particular case, IIj(cj) and ILj;(cj) are identical. This is a 
peculiar feature of the ID MI excitation spectrum, linked 
to the fact that in ID the hole and doublon cannot cross 
each other during their motion. Thus qualitative dif- 
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(d) 1D 



(b)2D 




(e)2D 



(c)3D 




(f)3D 
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FIG. 2. (color online). The imaginary part of the current 
[(a)-(c)], and the kinetic-energy [(d)-(f)] correlation functions 
in the bosonic Mott insulator for tn/U = 0.01 at zero tem- 
perature. The lattice constant and h have been taken to be 
unity. For ID [(a) and (d)], these two correlation functions 
are identical while they are qualitatively different for 2D [(b) 
and (e)] and for 3D [(c) and (f)]. 



ference must appear in the 2D and 3D cases. We thus 
compute IIj(cli) and ITjp(w) for 2D and for 3D by using 
an diagrammatic approach. We consider a doublon and 
holon with an infinite repulsive interaction which imple- 
ments the constraint that the two particles cannot be at 
the same point except when they recombine. For ID this 
method is in full agreement with (l(3 . [Hj]. The result 
is also shown in Fig. [2j One clearly sees the difference 
between the two modulations. Note that the anomalous 
structures in the amplitude modulations are related to 
the Van Hove singularities in the density of states. 

In summary we have proposed in this Letter to use 
small periodic phase modulation of an optical lattice to 
probe for the current autocorrelation function. The con- 
sequences of the modulation can be measured either by 
probing the absorbed energy of the system (for bosons) 
or by measuring the production rate of doubly occupied 
sites (for fermions). Such a phase modulation probe gives 
direct access to the frequency dependent conductivity of 
the system. 
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work was supported by the Swiss National Foundation 
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